^ 91-28105 I 

AEROSOL SPECKLE EFFECTS ON ATMOSPHERIC PULSED LIDAR 

BACKSCATTERED SIGNALS 

^ 32 ^ 



S.R. Murty 

Alabama A&M University 
School of Engineering and Technology 
Normal, Alabama 


The effects of refractive turbulence along the path on the aerosol speckle field propagation and 
on the decorrelation time are studied for coherent pulsed lidar systems. 


Introduction 

Lidar systems using atmospheric aerosols as targets exhibit return signal amplitude and power 
fluctuations which indicate speckle effects. 1,2 The speckle effects are manifested statistically as 
Rayleigh-distributed amplitude and exponentially distributed power of the backscattered signals. The 
signals from the various scatterers which make up the target combine to form a speckle pattern at any 
point along the backscattered beam. The speckle size is « A z/b, where z is the range from the target, A 
is the wavelength, and b is the radius of the illuminating beam on the diffuse target. Therefore, as the 
beam size on the diffuse target increases, the speckle size at the receiver decreases. Since the field across 
a speckle pattern is produced by the addition of a number of signals of random phase, the speckle pattern 
changes as each particle moves, thereby generating spatial and temporal changes in intensity. The 
particles move primarily due to atmospheric turbulence, and the changes in the speckle pattern are 
caused by the velocity turbulence in the illuminated region of the atmosphere which contributes to the 
backscattered signal at the receiver. For a pulsed laser, the backscatter region of the atmosphere is 
ctp/2 around the range point where t p is the pulse duration and c is the speed of light. 

In addition to producing signal intensity fluctuations, aerosol-generated speckle causes a 
decorrelation of the backscattered signals, as the aerosols are free to move about relative to each other 
along the beam axis and transverse to it. The resulting decorrelation time determines the temporal 
window of the return signal from which the power spectrum can be calculated to obtain wind speed. 

A laser beam experiences fluctuations in intensity and a reduction in coherence due to propagation 
in the refractive turbulence of the atmosphere between the telescope range point as shown in Fig. 1. 
Consequently the received signal characteristics are affected by a combination of refractive turbulence 
and speckle effects. The objective of this paper is to examine the effects of refractive turbulence along 
the path on the aerosol speckle field propagation and of the decorrelation time. 
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Chumside and Yura 3 have developed the basic theoretical formulation of the spatiotemporal 
correlation function for the aerosol return-signal field. Ancellet and Menzies 4 presented results of 
measurement using a pulsed TEA CO 2 laser and found a decorrelation time of 2-2.5 n s over a 1.3 -km 
range. Laser pulses shorter than ~ 150ps are broadened due to atmospheric turbulence, but the pulse 

spreading and wander in turbulence decrease as the pulse becomes longer. 5,6 In this work, the pulse is 
represented as a Fourier integral, and we make use of the two-frequency mutual coherence function 
to study propagation of the aerosol speckle field. 



Analysis 


Fig. 1. Sketch of pulsed lidar propagating through atmospheric 
turbulence using moving aerosols as a target. 


The transmitter and receiver are assumed to be colocated having a common telescope aperture 
diameter d. The laser transmitter is assumed to be pulsed with a pulse duration t p and a transmitted 
power Pt (t) at time t. The complex field of the pulse propagating in the turbulent medium at any range 
point z is given by 


Ui(z,v4) = J “ A(a>)uj (pjc,t) exp(i(ot)d a), (1) 

where A(co) is the Fourier spectral amplitude, uj is the monochromatic field amplitude at a distance z 
from the transmitter, p is a 2-D transverse vector in the aerosol target plane, and k is the optical wave 
number. Using the paraxial approximation and the extended Huygens-Fresnel principle, 7 u/ can be 
expressed as 

u/(pM = / d 2 ru T (r,t- Z -) 


x exp 



(P “ r) 2 + xp (r,pjc,t) . 


( 2 ) 


In Eq. (2), r is a transverse vector in the transmitter plane, y>(r,p Jcjt) is the sum of the logamplitude and 
phase perturbation experienced by a spherical wave propagating from r to p through the atmospheric 
turbulence, and uj represents the field amplitude at the transmitter given by 


Ut (r,t) = 

2 P t (/)] 
2 

Vz 

exp 


VfH 


71 a 



/ j 


where a is the beam radius at exp(-2) intensity, and f is the geometric focus of the telescope. 
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The complex field Ui (z,p,/) is assumed to be backscattered at range z by an atmospheric aerosol 
present in the illuminated volume. The scattered pulse propagates back to the receiver through 
atmospheric turbulence, and its field is given by 


U s (q,/) = f * A(o))u s (q Jc,t) exp(w)t)doj, 

V —oo 

where the monochromatic wave field u s is given by 9 
u s (q Jc,t) = “i (p&* -7) ) 


(4) 


x exp 


Ua + ^(q - p ) 2 +V'(p,q^) - 


(5) 


In Eq. (5), q is a transverse vector in the receiver plane and S(k) is the monochromatic backscatter 
coefficient of the aerosol. 


The pulse is assumed to have a center carrier angular frequency a>o and optical wavenumber ko 
and is narrowband. It is convenient to write Eq. (4) as 


U s (q,f) = exp[t(G><y - k 0 z)] I " d(oA(a) + o>o) 

J — 00 

xu s ( q , o) + wo, t) exp(icot). 


( 6 ) 


During a small time interval of <1 ms, the aerosol position is assumed to change due to a radial 
velocity v z , a transverse velocity V 7 - to a new position given by the following relations: 

z x = z[t + r) + z(t ) +v z r; y 

Pr =p(/ + T) =p(0 +Vr r. 

The backscattered field at time t + z, is given by 

U s (q,t + r) = exp[/a>o(t + T -z./ c )] I ® A (0 + <o 0 ) 

•/ — 00 

x u s (q, co + a>o, t + x) 

x exp[ia>(r + r -zv c )do), (7) 

and the monochromatic wave field at time and t + r is given by 
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u s (q ,o,t + x ) = — J* dr*UT |r,f + r - ^ j 

(q - p») 2 J 


ik , .2 , 

2F (Pr " 0 + 


+q f (r,p T ^ +6r) +q , (p r ,q^ +r). 

The cross-product of the received fields at time t and t + r is given by 

U s (qi,/) l £ (q 2 ,f + t) = f " f" doo*4 e (m) 

V — 00 v —oo 

X U, (qi,ft>i,t)/f e (<W 2 ) Us ($2,0)2, t + 0 

x exp [i (w\ — o ) 2 ) t ], 

where A e (<o) =/4 (<w + <wo)- 


( 8 ) 


(9) 


Equation ( 9 ) will be ensemble-averaged over the propagation path for refractive turbulence and 
over the illuminated volume for velocity turbulence to obtain the pulse spatiotemporal correlation 
function. This process is represented by double angular brackets, and the correlation function of the 
received optical fields is given by 


«U S (qi,t)f£ (q 2 ,* +r)»= f” do i f w do) 2 A e ((o{)Ae (o) 2 ) 

V — 00 J — 00 

xr(qi,q 2 ,a>i,t» 2 ) 

x exp[i (o i — 0 ) 2 ) /], (10) 

where r(qi,q2,<ui,<W2) = <«$ (qi,ft>i f)«? (q 2 ,co 2 , t +r) >>s a two-frequency mutual coherence func- 
tion. The determination of the effects of aerosol speckle, refractive turbulence, and pulse shape are 
now dependent on the evaluation of the function T. 

Ensemble Averaging 

Following Churnside and Yura 3 for the averaging over velocity turbulence, the aerosol velocities 
are assumed Gaussian-distributed within the inertial subrange, and we obtain the expected value of the 
mutual coherence function as 

r(qi,q 2 ,<ui,a> 2 ) = f°° dv z f°° d 2 y T f dzd 2 p 

J — 00 v — 00 V 

xu s (qi,<wi ,0 u* (q 2 ,e > 2 1 +r) 

xP(v z )P(V T )/V. (11) 


307 


In Eq. (11), V is the illuminated volume, and the probability density functions of the velocity 
components are given by 

p ( v *) = ^ 2 = exp [~(v z - V z f/2 o\ ], 

P(Yt) = — 3 exp[-(V r -V t ) 2 /o\ ], 

ZJtCff 

where V z , Vr and oj, of are the means and variances of the wind components along the z axis 
and transverse to it. 

The ensemble average over path turbulence is given by 

( u s (qi,fl>i,/) u* (q 2 ,<L> 2 ,t + r) ) = ex p^ 2 i(ki -ki)z 

Q*n 

-2ik 2 v z r] J' J' d 2 r\d 2 r$iT (rut 

xur (r 2 ,t + x exp^(p -n ) 2 
+^(qi -P ) 2 -§r(P +V 7 -T — r 2 ) 2 
-$(?2 -p - Vtz ) 2 ] 

xH(q h q2,rur2MM), ( 12 ) 

where H is the fourth-order mutual coherence function given by 

H =(exppP (p,ri^i,r) +W (p,qi^i^) + W* (p+V r r,r 2 ^ + r) 

+W*(p +V r r,q 2 ^ +t)>. (13) 

The change in the aerosol position has very little effect on the refractive turbulence encountered 
along the path for small r and will be neglected in evaluating the mutual coherence function, which can 
be written as 8 

H = exp [-(V 2 )(Di 2 -£>i3 +£>14 +D 2 3 +D 34 — D 24 +D 34 ) 

+2 C Xn +2 C x „l (14) 

In Eq. (14) Dq = D( | r, - rj |) is the two-frequency wave structure function of a point source 
given by 10 
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D{p) =2ip/p 0 f i -2A?2/Q 2 , 
po - (0.545fcoCn z) * 5 . 


Q~ 2 - 0.39 kodL^z, A 12 = |*i -k z \/ko, (15) 

where Lq is the outer scale of turbulence, c is the speed oflight, and C 2 is the turbulence structure 
constant assumed uniform along the path.Cn; —C x ( |r« — r y |) is the two-frequency covariance of the 
logamplitude of a point source given by 11 


Cxij 


0.\'$2ji 2 k\kj z C 2 



UU *Jo 


x[«|(r, —rj) (1 -01] s>n 


w 2 r(l -r)z 
2ki 


x sin 


u 2 t { 1 —t )z 

2kj 


(16) 


The integrations of Eq. (12) are lengthy and cannot be performed in closed form. We present the 
approximate results in this work for the simplified case of a collimated beam and negligible transverse 
velocity component and consider points along the beam axis. It is also assumed that A 12 « 1, and we 
make use of a quadratic approximation for the structure function. With these simplications, the 
integrations of Eq. (12) can be carried out, and the result is 


( u s (o,a>i,t) a* ( o,(o 2 , t 


+r)} - 


2S 2 F 2 

a 2 z\ 1 +F 2 +2a 2 /pl)jz 


x “r) ***■(' +T ~t) 

x exp J i(o> 1 — <02) Td — (g>i — a> 2 ) 2 /Q 2 


- (ki -k 2 )p 2 /2z -2 fe 


■)w 


1 + 


^k(/t 2 ^ 


\ 22 ) 


2ik\z -2ik 2 zr I-, 


} 


(17) 


where F = ko a 2 /2z is the Fresnel number of the source, Td si the time delay given by 


t ~JL 

d cko 


t 


2 /_? \ 


1 + CL / p§ 


1 + F 2 + 2a 2 /po 


(18) 
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the coherence bandwidth is given by 


«" 2 «0.39 Cltfz/c 2 , (19) 

and L 0 is the outer scale of the turbulence, and c is the speed of light. It may be noted that the 
pulse broadening is inversely proportional to the coherence bandwidth. We will now cariy/ out integra- 
tions indicated in Eq. (11) over the velocity spectrum and sensing volume V and replace S N/V by the 
aerosol backscatter coefficient /?, where N is the total number of aerosols within the volume V. The 
result is 


T(o,o,<o h o>2,) = | ~Pt (/|) Pt (' +r ] 14 

x exp [i (a>i —coi) Td -{<0 \ —a>i) /Q 

+2i(k\ -ki)z - Tiki V z a z x - 2o| k% r 2 ] dz. (20) 


The integrations over o\ and a>2 are carried out for a Gaussian pulse with frequency spectrum 
given by 


A(a>) =(2jioi ) 1/2 exp[(a> -a>o) 2 /2ai ], 

where Ou> is the bandwidth of the pulse. The variables are changed using 2zo'o — Q>\ + a>2, 
<5 = (o i — a >2 ■ 


The spatiotemporal correlation function is calculated for a Gaussian transmitted power distribu- 
tion given by 

= vi~ ex P 


where E is the total energy in the pulse, and the pulse width is related to bandwidth as 
t p =V2a < ^ 1 . The integration over z is performed designating arbitrarily the time of peak power by t = 
0, and the result is 
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«U S (o,t)lT s (o,t +r) » = 


cfiE Q T e 

2zo (Q 2 + 4oZ>)^ (Te + 


1 +6 F 1 +2 a 2 /p\ 


exp 


[-2<r|^r 2 


(* 2 +2Fz o z t/c ) 2 

v [ V , r 2 /4 + (7^ 4 - F z r/ c) 2 
|_4/ 2 r 2 

-i (tyox - + ^ovzr)] , 


( 22 ) 


where r? = oZ 2 + 4Q -2 and zo is the range. Equation (22) is the main result of this paper, which 
exhibits the effects of refractive path turbulence on pulse spreading and aerosol speckle. The first term 
in the exponential gives the decorrelation effect due to velocity turbulence along the beam of the axis; 
the second and third terms give the effect of refractive turbulence and velocity turbulence on pulse 
broadening and its contribution to the decorrelation of aerosol scattered signals. 


Discussion 

The spatiotemporal correlation function obtained by Churnside and Yura 3 has been extended to 
include the spectral effect of pulsed signals and refractive turbulence along the propagation path. An 
approximate result for a collimated beam along the beam axis is presented in Eq. (22) for the combined 
effects assuming a quadratic approximation for the structure functions. These approximations permit 
us to bring out explicitly the effect of time delay Td and coherence bandwidth Q on the correlation 
function. 

We consider a long path wind-measuring pulsed lidar having a pulse duration of 100 ns operating 
at 9.1 fx m wavelength. The path length is assumed as 100 km through a turbulent atmosphere with 
Cn = 6.4 x 10 -15 m~^ corresponding to a clear day in the lower atmosphere with an outer scale length 
Lq = 100m and o z = 1 wi The beam radius at exp(-2) intensity is taken as 1 m, which gives the Fresnel 
number of the source as F = 3.45. The transverse coherence length for this path ispo — 0.0465m. 

The coherence bandwidth Q and the time delay Td corresponding to these parameters are 
Q=4.1 xlO 11 rad/sand Tf = 1.65 x 10 -14 s.Using these values, we obtain o m =1.4xl0 10 rad/sand 
z e — 70.9 ns. These time scales are shorter than the correlation time of the order of l//s corresponding 
to moderate turbulence-induced aerosol dephasing, and the effect of pulse broadening is negligible for 
this case. When the correlation time due to aerosol dephasing approaches the nanosecond range, pulse 
broadening effects need to be considered. 

This work is supported by NASA grant NSG 8037. It is a pleasure to acknowledge the discussions 
with James W. Bilbro of NASA Marshall Space Flight Center. 
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